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This paper revisited quantum tunneling dynamics through a square double-barrier potential. 
We emphasized the similarity of tunneling dynamics through double-barrier and that of optical 
Fabry-Perot (FP) interferometer. Based on this similarity, we showed that the well-known resonant 
tunneling can also be interpreted as a result of matter multi-wave interference, analogous to that of 
FP interferometer. From this analogy, we also got an analytical finesse formula of double-barrier. 
Compared with that obtained numerically for a specific barrier configuration, we found that this 
formula works well for resonances at "deep tunneling region". Besides that, we also calculated 
standing wave spectrum inside the well of double barriers and phase time of double-barrier tunneling. 
" The wave number spectrums of standing wave and phase time show another points of view on 

resonance. From semi-numerical calculations, we interpreted the peak of phase time at resonance 
as resonance life time, which coincides at least in order of magnitude with that obtained from 
| uncertainty principle. Not to our surprise, phase time of double-barrier tunneling also saturates at 

long barrier length limit I — > oo as that of tunneling through a single barrier, and the limits are the 
■ same. 

O: 

in 



o 



PACS numbers: 



O-i I. INTRODUCTION 

Optical Fabry-Perot interferometer has a broad application in optics and has been already used in many optical 
^ | apparatus, such as optical filter, spectrometer, single-mode resonant laser cavity, etc.. The most central feature of FP 
O" 1 . interferometer is its excellent filter effect. For example, a commercial FP interferometer can easily achieve a finesse 
of about 200 (though may not be fabricated with a two parallel plates configuration, the principle beneath are the 
. same), and recently a finesse up to 130000 FP fiber interferometer was reported[l[. The high finesse of FP cavity 
' is accomplished through the multi-beam interference effect. On the other hand, interference is a crucial feature of 
quantum mechanics, so one may naturally wonder whether it is possible to construct a matter wave FP interferometer, 
r*"** ■ The answer is yes, and there have already been several mesoscopic devices with operating principle similar to that of 
FP interferometer available in the laboratory, thanks to the advanced mesoscopic fabrication technology 0- 

The foundational design responsible for the realization of high finesse FP interferometer we think is the resonant 
cavity formed by two parallel plates. To conceive a matter wave analogy, we think the well-known resonant tunneling 
is the most suited concept for this aim. Though we thought we had got an original idea of a matter wave FP 
interferometer, later we learned that Chamon et. aLhave already proposed such an idea much earlier, but with much 
more sophisticated concepts and interesting physics involved^- In their work, the tunneling particle is not non- 
interacting real particle, but low energy excitation-quasiparticle instead. And the double barriers we discuss below is 
replaced by two point contacts Q. Their proposal is later called electronic or Hall FP interferometer and has received 
much more attention^ recently. However, as our results is from a close analogy with optical FP interferometer, which 
require little calculation of Schrodingcr equation provided we have known transmission and reflection coefficients of 
single barrier in advance, and with this optical analogy, we have got a finesse formula matching well with the numerical 
results for resonances at "deep tunneling region" , we think these results may be still valuable. In addition to calculate 
parameters relevant to double- barrier, the matter wave analogy of an optical FP interferometer, we have also obtained 
the tunneling time of this barrier with phase time approach. We find that this time saturates at the same limit as 
that of single barrier, and we interpret this quantity at the local maximum as the lifetime of resonance. Compared 
with numerical results obtained from uncertainty principle, we think this interpretation is reasonable and may find 
application in the calculation of decay rate of resonances, e.g., in cosmology [2f|. 

The paper is organized as follows. First, we give a brief introduction to the concept of quantum tunneling and FP 
interferometer in section[Hj emphasized the similarity involved in between. Then we present a calculation of tunneling 
through a double-barrier following the standard quantum mechanical approach, i.e., solving Schrodinger equation. 
Next we calculate the transmission and reflection coefficients by a close analogy with that of optical FP interferometer. 



'Electronic address: blucscacat@126.com 



2 



(a) Square potential barrier 



(b) Double Square potential barriers 



FIG. 1: (a)Square potential barrier with heights A = 3.2eV, length / = 4nm; (b). Double square potential barriers with heights 
Ai = 3.6eV, length h = 2.5nm and A2 = 3.4eV, length h = 3.5nm, two barriers are separated by a distance d — 6nm. 



All these are given in section IIII1 In section IIV1 we calculate the phase time of double-barrier potential. With our 
calculation, we show that numerically, one can associate phase time at resonant wave number with the lifetime of 
resonance reasonably. At last, we give a conclusion in section IVl 



II. QUANTUM TUNNELING 



In this section, we give an illustrative example of quantum tunneling, which shows some similarity with that of 
optical FP interferometer and also inspired us to find a much resemblant matter wave counterpart. Consider the 

problem of a beam of particles (e.g., electrons or neutrons) with energy Hu)(k) = ^ scattered by the simplest 
square potential barrier of the form V(x) = A[Q(x — a)Q(b — x)] (see Fig 1(a) ), 
the solution of the corresponding Schrodingcr equation is 

$(x,t) = <j)(x)e-™ {k)t (1) 

fh{e lkx +Kc-' lkx ), x<a; 
with (j){x) = { ^J%{ce fjx + de-^), a<x<b; (2) 



%Te lkx , x > b. 



Where /3 2 = j?r(A — fiw(fc)), I is the incident flux, and the transmission and reflective coefficients are given by 

[(f - f ) sinh[/3(6 - a)} + cosh[/?(6 - a)} 
4(| + |)sinh[/?(6-a)]e 2lfe 



T= —p >,., rn „ — ^> (3) 



n = — si 'i (4) 

i(f - f ) sinh[/3(& - a)} + cosh[/?(6 - a)} 
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It is interesting to note that if we make a sign change, i.e., A — > —A and /3 — > i/3, then potential barrier becomes 
potential well. Using the well-known relations 

sinh(i/3x) = isin/3a;, cosh(i/3x) = cos/3x, (5) 

the transmission rate (the square modulus of T) now becomes 

2kPe~ lkL l2 1 



1 2fc/?cos[/3L] — i(/3 2 + k 2 ) sin[/3L] 1 1 + I (| - 2)2 sin 2 [/3L] 



with L — b — a. (6) 



We note that the transmission rate is a periodic function of j3L. Under certain circumstances, the transmission rate is 
equal to 1, in other words, the incident beam of particles completely transmit through the well. This is the well-known 
resonant transmission phenomena [f|. 

This phenomena bears some similarity with a photon beam transmitted through a Fabry-Perot (FP) interferometer, 
though in the latter case the transmission rate is a periodic function of photon energy E\ or wave number fc, i.e., 



(7) 
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with 

3 = tiV#/(1 - 91), E x = hu = hc/X. (8) 

^ defined in ^ is the finesse of FP interferometer, and only depends on reflectivity in this case. Note that if 
1 — 3$ <C 1, the transmission rate is very small except that incoming photon is at resonant energy E\ = ch ^ (where 
m g Z); while for a quantum well, it was supposed to be 1 in classical logic but is shown to be less than 1 in 
quantum mechanics in general except at resonances. Of course, as a quantum analogy of optical FP interferometer, 
quantum well indeed captures the essential concept (resonant transmission) and many interesting works have devoted 
to this problem^, however, we do not think this analogy is the most suited one. We find maybe tunneling through 
a double-barrier is a more suited analogy for this study. In the next section, we will calculate the tunneling rate 
through a double-barrier configuration, in which the physic picture bears much resemblance than that of a quantum 
well to the optical FP interferometer's. In addition to the standard quantum mechanical calculation, we give another 
derivation of the transmission and reflection coefficients based on the close analogy between optical FP interferometer 
and a double-barrier potential. 



III. DOUBLE-BARRIER POTENTIAL 



A. General Result 



The solution of a beam of particles scattered by a double-barrier potential 
U db (x) = Ai[9(x-ai)6(6i - x)] + A 2 [Q{x - a 2 )@{b 2 - x)], with b { 
can be obtained by solving the corresponding Schrodingcr equation 



dt 



~^+u db {xm x ,t), 



li, i = 1, 2; a 2 — b\ = d. (9) 



(10) 



where d is the distance between double potential barriers and l\, l 2 and A±, A 2 are the widths and heights of 
the corresponding barriers respectively. The double-barrier configuration is shown in Fig 1(b) The two-barrier is 
reminiscent of the parallel two plates in FP interferometer, later we will use this resemblance to derive the transmission 
and reflection coefficient in another way. 
With the given ansatz 



^{x,t) = (!>{x)e- lu > Wt 

' z lkx + K db e 



with 



Ci& 



(*) = 



ikx 

die-? 1 *, 

ikx 



ae + 7e 
c 2 e hx + d 2 e~^ x , 
T db e ikx , 



x < a\\ 
ai < x < b\; 
b\ < x < a 2 ; 
a 2 < x < b 2 ; 

x > b 2 , 



(where we have normalized incident flux Iq to 1) and the continuity condition of wave-function 



<p(a, - 0) = <f>(ai + Q), 



d(j>(x) 



d(f)(x) 

\x— ai— = ] 

dx dx 
we can obtain the transmission coefficient 

e -ik(d+h+l 2 ) 



\x— ai+0: 



0) - 0(&, 
dx 



0), 



dcf)(x) 
dx 



|x=a<+0, ( i = 1 > 2 ) 



77, 



db 



F[k-p u h-d-hMY 



[cosh(/3 2 /2) 



e lM [Mi sinh(^iZi)][M 2 sinh(/3 2 Z 2 )] + e- ifc(1 [cosh(/3iZi) + iK x sinh^iZi) 
- iK 2 svah(p 2 l 2 )], 



where M,;, Ki in (pT6|) are defined as 

Mi = 



2 l k 



+ ■5"). K * 



2 { k 



-); (i = l,2). 



(11) 
(12) 



(13) 
(14) 

(15) 
(16) 

(17) 
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FIG. 2: A matter multi-wave interference for a beam of particles tunneling through a double-barrier potential. 



The other coefficients are given explicitly in Appendix A. 

We have already observed the formal analogy between two barriers with the two parallel plates in FP interferometer. 
Next we need to know the reflection and transmission rate and the corresponding phase changes for a particle tunneling 
through or reflected by each barrier respectively. These have already been obtained in section UH here we decompose 
the amplitudes 

T * = \n = TT-TTo ■ 1 1 2„,, 11 /^ (18) 



M?sinh 2 CSA)]V2' 
MiSmh(Pik) 
[1 + M?smh 2 ((3 t h)} 1/2 



R * = I'M = M , „.r9 ■ ; 2,n i vTTTo ' (* = 1' 2 ) ( 19 ) 



and phases 



ki = -kk - <j>i, (20) 
b ri = -it/2 - fa, fa = arctan[A, tanh(ftL)]; (* = 1,2) (21) 



of the corresponding formulas (|3l4p explicitly. Next we will use this optical FP analogy, or more precisely multi-wave 
interference (shown in Figj2j) to regain the transmission and reflection coefficients of double barriers. 

For a beam of particle scattering on a double-barrier potential, it will be partially reflected and partially transmitted 
through the first barrier. The reflected component is denoted as n; the transmitted component propagates forward 
through a distance d and meets the second barrier. Then it will partially transmitted and partially reflected again, 
the transmitted component is denoted as t\, while the reflected component propagates backward through another 
distance d and then meets the first barrier. Then it will partially transmitted with transmitted component denoted 
as r 2 , and the reflected component will repeat the story again and again, as shown in Fig(21 The final amplitude of 
transmitted and reflected waves are the summation of these transmitted and reflected partial waves, i.e. 

T db = Y^k = TiT 2 exp[i(&i + & a ) + ifc(d + h + h)} + T X R 2 R{F 2 exp[i(<f> tl + <f> rl + fa 2 + fa 2 ) + ik(3d + h + l 2 )] 

i=l 

+ Ti(R 2 Ri) 2 T 2 exp{# tl + 2(fai + fa 2 ) + fa 2 ] + ik(5d + h+ l 2 )} + ... 
= T X T 2 exp[i((f> tl + <j) t2 ) + ik(d + h + l 2 )] 



1 - RiR 2 cxp[z(0 r i + fa 2 + 2kd)} 

CXp{z[0 t i + fa 2 - ((f>rl + fa-2) + Hh + h)}} 



^} j? 2 cxp(ikd) — jt^t- exp[— i{4> r i + fa-2 + kd)] 
From (|18I19|) and (|20I21|) , we can obtain 

= A/ismMftiOA/asinh^iia), - ^H^' + ^I = {[i + M \ sinn 2 ^)] [1 + Af| sinh 2 ^)]} 1 / 2 exp[t(<£i + <h)] (23) 
and the phase factor 

fal + <Pt2 ~ (fal + 0r 2 ) + Hh + h) = 7T. (24) 

Substituting (|23|24j) into ([22} and utilizing the identity 

exp(ifa) = [1 + Mf sinh 2 (/3J l )] 1/2 cxp{i arctan[A', tanh(^/,)]} = cosh[/3 4 L] + iK l sinh[/3;L], (i = 1, 2) (25) 

J- i 



(22) 
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we can obtain the transmission coefficient 

Tdb = y^Kh^hM- (26) 

Compared with (fT5|) . there missed a phase factor e ~i k i d + l i+h) _ This phase factor can be inserted back from the 
comparison with single barrier case ([3j in the limit d — > (which picks up a phase factor e _lfc ( ( i+'2)) an( j that with 
right moving plane wave e lkx in the thin barrier limit l\ — I2 — > 0. 

Similarly, the reflective coefficient can be obtained with the same procedure utilizing multi-wave interference, and 
the result coincides exactly with (|63[) calculated from Schrodinger equation. For details, see Appendix B. 

From the derivation of transmission coefficient (|2l)]) . we can immediately get the resonant tunneling 5] 6] [20] condi- 
tion by observing that 



1 |2 (1-JZiiZa)- 



1 - RiR 2 exp[i(q> rl + <f> r2 + 2kd)] 1 + (1 t\% Y 2 sin 2 [fcd + ±(</> rl + r2 )] 
hence the transmission rate gets the maximum value 



(27) 



T m ax = {- 5 5 ) (28) 

1 — r£irt2 

only when 

kd + — ((f> r i + 4>r2) = mir, m G Z. (29) 

Note this condition, i.e., $(fc) = kd — \{<Pi + 02 + t) = m7r, m £ 2 is obtained by a direct analogy with the 
expression of FP interferometer, ([7]). From the same formal analogy, we can also define the finesse of a double-barrier 
interferometer as 

= n^Rjh _ 7r[A/ 1 sinh(/? 1 ^ 1 )M 2 sinh(^ 2 / 2 )]^{[l + M 1 2 sinh 2 (/? 1 ^ 1 )][l + M 2 2 sinh 2 (^ 2 / 2 )]} 1 / 4 
dh I-R1R2 {[l + M 2 sinh 2 (/3 1 ; 1 )][l + M 2 sinh 2 (^ 2 ; 2 )]}V2_M lS i n h( / 3 1 ; 1 )M 2 sinh(/3 2 Z 2 ) ' 

This formula works well especially for resonance at "deep tunneling region" , i.e.. the peak of resonance in wave number 
spectrum is very sharp. This will be confirmed numerically in below. From (|28l) we know in general T max < 1, except 
Ti = T 2 . The simplest case where this condition is fulfilled is the symmetric double-barrier. We will utilize this 
special case (i.e., l\ = I2 = I and A\ = A 2 = A) to give a numerical calculation of finesse and resonance condition in 
the following. The numerical calculation in turn serves as a cross check of the usefulness of our analogy. 

B. Symmetric Double-Barrier 

In this case, the transmission coefficient is 

e -ik(d+2l) 

Tdbs = F[k,p,l,d]' (31) 

with 

F[k, /3, 1, d] = [M sinh(/3Z)] V M + [iK sinh(^) + cosh(/3/)] 2 e~ ife<i , (32) 

where 



The transmission rate is given by 



Tdbs ~ lTdbsl2 -\T[k,p,l,d}\2> (34) 
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FIG. 3: left figure shows transmission rate with respect to the incoming particle wave-number k. The green curve is plotted with symmetric 
double barriers: heights A — 10.36cV, length / — 1.2A and distance d — 6.5A. The blue one are two identical curves plotted with pair of parameters 
{A\ — 10.36eV, li — 1.2A) and {A2 — 9.6cV, I2 — 0.6A) interchanged, the distance between two barriers is also d — 6.5A. The middle and right 
figures with cross points corresponding to the roots of two transcendental equations of $391 and J40B respectively, and those indicated by dark disks 
are roots of the simultaneous equations. Note horizontal axis of wave number is plotted with unit k e — \/2m e E s /h, i.e., the wave number of an 
electron with kinematic energy E s — lcV. 



with 



\T[k,/3,l,d]\ 2 = 1 + M 2 [cosh(2/?0 - 1] 



M 2 



M 2 M 2 
1 + — [cosh(2/3Z) - 1] + A' sinh(2/3Z) sin(2fcd) + {cosh(2/3i) — [cosh(2/3i) - l]}cos(2fcd) 



1 + M 2 [cosh(2^) - 1] { 1 + — [cosh(2/?0 - 1] \ [1 + sin(2fcd + 5)}, 



cosh(2/3Z) - M_[cosh(2/80 - 1] 
j£Tsinh(2/3Z) 



where in the last step of ([35)1 we have defined a fc-dependent angle 

tan S 

and utilized an identity 

M 2 

[A'sinh(2£Z)] 2 + {cosh(2^«) — [cosh(2/3Z) - l]} 2 = 



M 2 

1 + — [cosh(2/30 - 1] 



(35) 



(36) 



(37) 



From (|34|35p . it is easy to find that when 2kd + S = 2Nir + T dbs = 1, i.e., incoming particles completely 
transmitted. While for the case 2kd + S = 2Nir + |, the transmission rate is locally minimal, and corresponds to the 
minimal tunneling 



( 



I + R1R2' 



(38) 



in section UlIAl 

The above two local extremal conditions can be summarized as 



tan[<5(fc)] = cot(2fcd). (39) 
To determine at which wave-vector k there is a resonance, we have to solve simultaneous equations Q39p and 

Ksmh(2pi) 



sinf2fcdl = — - 



1 



Ml 
2 



cosh(2/3Z) - 1] ' 



(40) 



However, these equations (including the simultaneous equations to determine local minimal of transmission) are 
transcendental equations, so we can only solve them numerically. 

The results are illustrated in Fig[3l In the left subgraph of FigEl we have plotted the curves of transmission rate 
(I15p vs incoming wave-number k. The blue curve shows the transmission spectrum of general asymmetric barriers 
(with parameters (3\ 7^ /3a and l\ 7^ £2)- Note that this is actually a layered curve of two identical ones, with 
pairs of parameters A\ <-> A2, l\ <-> I2 interchanged simultaneously. This permutation symmetry is manifested in the 
transmission formula (|16p . and inherits from the left-right symmetry of 1-Dim Schrodinger equation (|10p . i.e., particles 
incoming from left must have the same transmission rate with those coming from right. As we already stressed in the 
end of section UlI A[ the maximum transmission rate in the asymmetric-barrier case is less than 1 in general. While 
for the symmetric case, the green curve, its peaks approach unit and provides a much ideal example of optical FP 
analogy. The two pairs of curves in the last two subgraphs of Fig 13] correspond to left and right hand side of equation 
(|39p and (|40p respectively, and the coincident cross points (marked with black dots) at the same wave number k in 
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Transmission -Rate Probability 




Wave - Number 



(a) Transmission Spectrum for Optical FP (b) Transmission Spectrum for Double Barrier 

Interferometer 

FIG. 4: Transmission rate plotted as a function of wave number k. (a). The relevant parameters for optical Fabry-Perot 
interferometer are 7? = 80%, do = 2cm and wave number k is in unit 10 6 /m. (b). The relevant parameters for a symmetric 
double-barrier are Ax — A2 — 10.36eV, li = h = 1.2A and m = m e , d = 7 A. Horizontal axis of wave number is also plotted 
with unit fc e = y/2m e E 3 /h. 



these two subgraphs denote the solutions of simultaneous equations, (|39I40[) . Thus these points also pin down where 
the resonances occur. We can utilize this pictorial method to obtain the numerical solutions of these transcendental 
equations, and analyze the finesse of each resonances (or filter effect) for a specific double-barrier numerically. 
The finesse of an interferometer at certain wave number is defined as 

3or = Ad f ree / AcTFWHM, (41) 

where Aaf ree is the distance of nearby resonances in the wave-number axis, and AapwHM is the full width at 
half maximum (FWHM) of a resonance (peak in the fc-spectrum). For convenience, we plot transmission spectrum 
of optical FP interferometer and symmetric double-barrier in FigfJ] as a direct comparison. For a R = 80% FP 



interferometer shown in Fijj|4(a) the finesse defined in (|8|) can be easily worked out as $ = 14.0496. While for a 



double-barrier matter wave interferometer, we have to work out the wave numbers at which Tdbs = \ is satisfied for 
each resonance, and this is determined by the transcendental equation 

M 2 [cosh(2/?0 -1] ll+ — [cosh(2j8i) - 1] \ [1 + sin(2fcd + 6)] = 1. (42) 



Wave number at resonance 
Aa di3 


0.742007 
0.742007 


1.47909 
0.737083 


2.20664 
0.72755 


2.92188 
0.71524 


0.735410 
0.735410 


1.467477 
0.732066 


2.19272 
0.725246 


2.907172 
0.714449 


2 at left side 
i at right side 
Act fwhm 


0.739935 
0.744104 
0.004169 


1.4698 
1.48866 
0.01886 


2.18166 
2.23328 
0.05162 


2.86621 
2.98555 
0.11934 


0.731927 
0.738966 
0.007039 


1.450947 
1.484995 
0.034048 


2.145098 
2.247656 
0.102558 


2.793121 
3.077995 
0.284874 


ddb 


175.215 
177.562 


38.7312 
38.8839 


14.1509 
14.0283 


6.12092 
6.05179 


103.252 
104.338 


21.3461 
24.684 


7.08665 
9.99817 


2.55128 
4.86271 



TABLE I: This table is divided into two parts. The left four columns are the data corresponding to symmetric double barriers 
shown in Fig 4(b) while the right four columns are data corresponding to asymmetric double barriers shown in Fig 5(b)| The 
first line in this table are wave number at which resonances occur. The second and third lines are wave numbers at which half 
maximum is approached from left and right in fc-axis. The fourth line are FWHM corresponding to resonances showing in the 
first line. The last two lines are finesse obtained with numerical method and the analytical expression (|30[) respectively. All 
these numbers are in unit k e = V2m e E s /h « 5.12289 x 10 9 /m. 



For a specific double-barrier shown in Fig |4(b)| we work out the wave-numbers k where the resonances occur and 
also the roots of transcendental equation (|42|) using the pictorial method shown above. The results are summarized in 
Table[l] Since there are only four peaks in tunneling region, we only need to give 4-column of results. The approximate 
linearity of the phase $(fc) = kd — \ { <j>\ + (f>2+ tt) factor in (|2"7|) shows that the occurrence of resonances is nearly 
periodic with respect to k (see Fig 5(a) ), thus Aa,i is (k) varies slowly and we can approximate the free distance between 
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TransmissionRate 





(a) Phase function <E>(fc) 



1.0 1.5 2.0 2.5 3.0 3.5 

(b) Spectrum of Asymmetric Double-Barrier 



FIG. 5: (a). An approximate linear curve of $(fc). The five dashed horizontal lines correspond to 0, n, 2ir, 3ir, 4ir, while the 
cross points of these lines with the blue curve show out at which fc, the resonant tunneling condition (|29|) is satisfied. The other 
parameters specify the double-barrier is the same as in Fig |4(b)[ (b). The transmission rate for an asymmetric barriers, with 
Ai = f0.6eV, h — 1.5 A; A 2 = 8.7eV, h = l.A and m = m e , d — 7 A. It also show nearly four peaks in the tunneling region. 



nearby resonances with their averages, i.e., 

I 4 

A<7 free = T ^ A(X dis = 0.73047. (43) 

i=l 

From table U we see that Actfwhm increases with increasing wave number, hence $ or decreases with fc increase. In 
the last two rows, we see that the finesse obtained from analytical expression (|30[) . ^db, matches well with $ or obtained 
numerically for the symmetric double-barrier. Even for an asymmetric double-barrier, we see the results match well 
for resonances in "deep tunneling region", i.e., the fifth and sixth columns of data in table HI 



IV. RESONANCE AND TUNNELING TIME 



Before investigating resonance behavior, first we consider the standing wave phenomena happening inside a FP 
cavity and inside a well formed between two barriers. We will sec this phenomena has close relationship to resonance. 
For a FP interferometer, standing wave forms due to the superposition of left moving and right moving multi-waves. 
The relative intensity is given by 



l + &-2V&cos(2k(d -x)) rr 
I[k,x) = — — Tpp, (44) 

where the origin of position x-axis is assumed to be at the edge of left mirror of FP cavity. I(k, x) as a function of k 



is shown in Fig 6(a) We see that except certain missing peaks, the peaks appear nearly at the same position fc — 



do 



(m G Z) in fc-axis. This fact is due to that I(k,x) is a product of Tpp and a pre-factor, and the latter in general 
only causes a slight shift of the peaks in fc-axis. The consequence of significant modifications caused by the pre-factor 
are the missing peaks. These missing peaks are due to the presence of wave nodes at certain fc, and reflect most the 



standing wave nature in the cavity. A more complete graph of I(k,x) is shown in Fig 7(a) in Appendix B. 
Similarly, the relative probability of finding a particle in the free region (x € (61,02)) between two barriers is 

Pb 1 <x<a 2 (k,x) = \4>{x)\ 2 

= {1 + A/|[cosh(2/3 2 / 2 ) - I] + M 2 [sinh(2/3 2 /2) sin[2fc(a 2 - x)] - J<- 2 [cosh(2/3 2 / 2 ) - I] cos[2fc(a 2 - x)]]}T db . (45) 

This time the standing matter wave is formed by the interference of left moving and right moving components, ae ikx 
and rye~ lkx . A 3-dimcnsional plot of Pb 1 <x<a 2 (k, x) is shown in Fig |7(b)| two slits of which arc given in Fig |6(b)| 
Similar to /(fc, x) in (|44j) . < x <a 2 {k, x) is also a product of transmission rate Tdb with a pre-factor, so except certain 
missing peaks, the relative probability approaches maximum nearly at the same fc where resonant tunneling happens. 
As is transparent in Fig |6(b)| the pre-factor in (|45|) is also responsible for the occurrence of wave-nodes (or missing 
peaks) and reflects the essential feature of standing wave. 

Note that the wave function in (jlip is unnormalized, thus there is no need for Pi, 1<x<a2 (k,x) to be less than 1, 
similar to the relative intensity of a standing wave in FP cavity. Further more, due to the constructive interference 
effect of the matter wave components, ae lkx and je~ lkx , Pb 1 <cx<a 2 (k,x) is much larger than 1 at certain resonant 
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4.0645 4.0646 4.0647 4.0648 4.0& 

(a) Optical FP cavity 



Wave - Number 



A 



Wave — Number 



(b) Double barrier potential 



FIG. 6: fc-spectrum of standing waves. The missing peaks are due to the presence of standing wave nodes, (a). The parameters 
for the spectrum of standing light wave inside a FP cavity are R = 80%, do = 2cm and k is plotted with unit 10 6 /m. The blue 
and red curves correspond to I(k,x) with position x = 5.03mm and x = 1cm respectively, (b). The parameters for the spectrum 
of standing matter wave of a double-barrier are A\ = 9.6eV, h = 1.2A; A2 = 25.8eV, h = 0.8A and m = m e , d = 7A. The 
blue and red curves corresponds to Pb 1 <x<a 2 (k,x) with position x = 3.25A and x = 4.6 A respectively. Wave number is in unit 
fc e as mentioned above. 



wave numbers. On the other hand, there is a flux conservation condition 

ih 



2m' 



(x)'(j)(xy - (j>{x)*' (t>{x)] = Const, 



which insures that 



\% b \ 2 = |a |2 '- 12 



171 



\T db \ 2 + \K db \ 2 = 1. 



(46) 



(47) 



So the net flux in this 1-dimensional problem is a right moving one, together with the reflected left moving flux, the 
total flux is always conserved. 



Next we turn to the tunneling time of double-barrier, this may shed some new light on resonance effect. As a physical 
process, it must take certain time for a particle to tunnel through a barrier. However, time in quantum mechanics is 
a rather controversial subjects [Hj]. Up to now, there does not exist a well accepted definition on tunneling time in 
the literature, though there does have several well-known proposals [13 EH O and a tremendous amount of works on 



this subject [l5 1 16j 19j](for details, see [H) [14]). Here we do not try to give a self-contained discussion or clarification 



on this obscure subject, instead, we just try to utilize the phase time approach to estimate how much time a resonant 
particle spends inside a specific double-barrier. 

To make life simpler, consider the symmetric double-barrier, then the transmission coefficient pip can be rewritten 

as 



T d 



\Tdbs\exp[iS p (k;0)], 



(48) 



where S p = 9(k; O) — k(2l + d) and 6(k; O) is the phase shift for a particle tunneling through double barriers. Here 
O denotes all the other parameters (like barrier length I, distance between two barriers d) except the wave number. 
The phase shift can be written explicitly as 



9{k; O) = arctan[ 



[M 2 - K 2 cosh(2/3J)] sin(fcd) - Ksmh^fil) cos(fcd) 



cosh(2/3/) cos(fcd) + K sinh(2/3Z) sin(fcd) 



from which one can obtain the phase time 



_ 86(k;0) dk _ m d0(k;O) 
Tdbs ~ d{k) 'dE ~ lik dk ' 



(49) 



(50) 



The explicit analytical expression of r dbs is presented in Appendix C. Before analyzing r dbs , we think it is valuable to 
review the result of single barrier first, and then compare it with that of double-barrier. The phase shift for a single 
barrier is 



(51) 
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from which we can get the corresponding phase time 

m M 2 sinh(2/3Z) 



Tsb 



K{kl) 



Hk (3[1 + M 2 smh((3iy 



(52) 



We show in Fig 
I and barrier heig 



a) and Fig |8(b)] that how r s b varies with respect to its arguments, wave number k, barrier length 
it A. At much small length I, T s b decreases monotonously with increasing wave number. This is 



anticipated, since energetic particle moves faster than less energetic one, obviously it will take less time for an energetic 
particle to cross over the barrier. At larger length I, T s b decreases with increasing k for k < fc , and increases slowly 



for k > ko, where fco is a stationary wave number satisfying 



dr 3l 
dk 



k=k 



0. The barrier height dependence of T s b is 



quite similar. At rather small I, to cross over a higher barrier, a particle with a certain k will spend more time, while 
for larger I, phase time decreases with increasing A. Generally speaking, thin barriers slow down particle's velocity, 
while for thick barriers, the more repulsive the barrier is (i.e., the larger barrier area At), the less time it will spend 
in tunneling. In other words, t s & approaches ^(2 + A/E) when (31 <C 1 and saturates at the value ^ when (31 3> 1. 



The saturation behavior of T s b with length I can be apparently seen in Fig |8(a)| and is a signal of the well-known 
superluminal phenomena discussed in the liter ature 0, |sj|. However, casuality isn't violated as shown in[I3|[I3l- We 
suspend our review here and return to the discussion of double-barrier. 

At much small barrier distance d and small barrier area (i.e., small A and I), the phase time behavior of double- 
barrier with barrier length I is qualitatively the same with that of a single barrier with length 21. This similarity is 
not surprising since symmetric double barriers with barrier length / becomes single barrier with length 21 in the limit 
d — > 0. Continuity guarantees this resemblance in the case of small barrier distance d and barrier area Al. This can 
be confirmed from the comparison of Fig [9(b)] with Fig[~ 



This resemblance breaks down with larger barrier distance. In this case, the presence of a well inside two barriers 
renders this time behavior quite different from that of a single barrier. From the analytic expression 



1~dbs 



mV[k,(3,l,d] 



T, 



dbs ■> 



(53) 



we see the effective length of tunneling region is a product of tunneling rate Tdbs and T>[k, (3, I, d] (T>[k, /3,l,d] is 
defined in Appendix C), both of which contain periodic sine functions of kd, thus it is obvious that Tdbs must have 
certain resonant behaviors. This is visualized in Fig 9(a) with Fig |9(b)| The qualitative behavior can be summarized 
as T dbs ~ {d + 21 + 2ZM[| + M(kl)s'm(2kd)}}/v when pi < 1 and r dbs w 2/2 when (31 > 1 (where v = Hk/m). 
Interestingly, the phase time of both double-barrier and single barrier saturate at the same limit when (31 1, so the 
superluminal problem haunts us again in the double-barrier case. This is not surprising, since in the limit I — > oo, 
the presence of a finite well is no longer important, Tdbs thus approaches the same limit as that of T s b- We will not 
hesitate here about the subtle causality problem in quantum tunneling. Rather, we turn to the more practical and 
interesting analysis of resonant behaviors of Tdbs instead. 

As is apparent in Fig |9(a) the presence of a well is responsible for the landscape of peaks and valleys in phase 
time of double-barrier. Generally speaking, for a given double-barrier heights, a wider well allows the formation of 
more peaks, hence resonance. This fact is not only reflected in tunneling rate Tdb (or Tdbs) and Pb 1 <x<a 2 (k, x), but 
also in phase time ((53)) . So we guess the sharp peaks in tunneling time can be attributed to the longer life-time of 
resonances. As a numerical check, we use uncertainty principle 



St ~ H/SE 



(54) 



to give an estimate of the life-time of resonances for a specific double-barrier, since we have already calculated 
AaFWHM in table HI it is easy to obtain the width of resonance by SE ~ h 2 kSk/m, where we have approximated Sk 
with AaFWHM ■ The results together with that calculated from ([53]) is presented in Table HU 



Wave number at resonance 


0.742007 


1.47909 


2.20664 


2.92188 


0.735410 


1.467477 


2.19272 


2.907172 


AaFWHM 


0.004169 


0.01886 


0.05162 


0.11934 


0.007039 


0.034048 


0.102558 


0.284874 


2t uc 

Tdb 


213.041 
213.311 


23.6247 
23.8112 


5.78567 
5.95756 


1.88997 
2.08377 


127.31 
134.175 


13.1899 
15.1924 


2.93055 
4.33754 


7.95755 * 10 _1 
1.73741 



TABLE II: t uc is the life time calculated from uncertainty principle, while Tab is that obtained directly from phase time. As in 
the table HI the first four column are time quantities for resonance of symmetric double barriers while the last four column are 
that of asymmetric double barriers. The time unit is in femtosecond (10 _15 s). 
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We see that these two estimates match very well with an error less than 0.3fs for symmetric double barriers. Using 
semi-analytical method, we also obtain the phase time of an asymmetric double-barrier with parameters specified 



in Fig 5(b) We see that for high finesse resonance, the life time estimated by phase time approach and that with 
uncertainty principle coincide well at least in the order of magnitude. Of course, if the estimation of life-time is 
replaced by Tdb(k, d, h, l 2 ) — [T s b(k,h) + T s &(fc, 1%)], the results will be more close to 2t uc . The factor of 2 in 2t uc is 
inserted to fit the data well. Up to now, at least in the order of magnitude estimation, we do provide an alternative 
method, though we have not worked out why the phase time approach can be a good evaluation of the lifetimes of 
resonances in double-barrier. Probably the fact that the energy eigenvalue of bound states in a well with width d 
formed by two infinite thick barriers with height A satisfy M 2 sin 2 [kd] = 1 may provide a hint of the underline reason. 
The detail analysis of this problem is beyond the scope of this paper. 



V. CONCLUSION 



We have revisited the problem of quantum tunneling through a double-barrier with a close analogy of optical Fabry 
Perot (FP) interferometer. Though there have many interesting schemes of matter wave FP interferometer [U [HQ 21 



and the physics involved are much richer than double barriers, our derivation of transmission and reflection coefficients 
are based on a common concept between the two analogs, multi-wave interference. Not only our derivation relies only 
on a knowledge of reflection and transmission coefficients through two single barriers (without solving the Schrodingcr 
equation of double barriers), but utilizing the same analogy, we also obtain a finesse formula (|30p . which is in good 
accordance with the definition, (|41[) . for resonances in "deep tunneling" region, especially for the symmetric double 
barriers. 

Then we discuss the relevant standing wave phenomena, which is also a common feature of FP interferometer and 
double-barrier tunneling. Moreover, it is associated to the occurrence of resonance. Utilizing phase time approach, we 
calculate the phase time of double barriers in comparison with that of a single barrier. The supcrluminal phenomena 
or Hartman effect 0] also resides in the double barrier case, which may be a general feature of quantum tunneling. Of 
course, several s ugg estions have been proposed to mitigate or even eliminate this effect and the associated causality 
problems [3 Hal • Here we don't discuss this topic in this paper, we feel that a suitable solution should be applicable 
to various barrier shapes and may even involve relativistic quantum formalism. Tunneling time of various shapes of 
double barriers has already been discussed in the litcrature[23] [24] 25] and this paper only provides a slight touch on 
phase time of double square barriers. We find phase time of both symmetric double-barrier and single barrier saturate 
at the same limit when pi — > 00. Further we focus our attention on the peaks and valleys in the configuration 
of double-barrier phase time. With a numerical support, we interpret these peaks as results of longer life time of 
resonances formed inside the well between two barriers. The numerical results match well with that estimated from 
uncertainty principle, especially for resonance with high finesse. Thus we think phase time may provide a helpful 
method to the evaluation of resonance's life time. Furthermore, we even guess that other time quantities, like dwell 
time, Larmor time, etc. [Tl| [HI, may have interesting applications in the estimate of life-time of resonances in quantum 
tunneling. Though the underline relation of phase time and the life time of resonance is not clear, we think it is valuable 
and interesting to explore it in the near future. 
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VI. APPENDIX A 

Various parameters characterizing tunneling through double potential barriers 

U d b(x) = Ai[Q(x - a!)e(h -x)] + A 2 [e(x - a 2 )0(b 2 - x)], with h - en = h, i = 1, 2; a 2 - h = d, (55) 

with wave number k = ^ 2 ™ E and = V 2m (^' j§ ^ _ ^2) is given in this section. To make the expressions as 
compact as possible, we define 
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Then the transmission coefficient Tdb in (fTS"j) is 

e -tk(d+h+l 2 ) 

Tdb = TfcpufcAhMY (57) 

F[k- p u fo; d; h,h] = e ikd [Mi sinh(piZi)][M 2 sinh(/3 2 ? 2 )] + e.-' iM [cosh(/3iZi) + iK x sinh(A/i)] 

[cosh(/3 2 Z 2 ) + iK 2 sinh(/3 2 / 2 )]. (58) 

The other coefficients in (fT2"]) are 

C 2 = ! e -fe(«Wl+W(l + ^)^-l e (*-ft)^, d 2 = I e fe(^i+^)(1 - ^)^-l e (*+/3 2 )a 1; (5g) 

2 /? 2 2 /3 2 

a = e- lfc ( d+il )[cosh(/3 2 Z 2 ) + iK 2 smh(0 2 l 2 )\ JF" 1 , 7 = -ie lfe ( d + il+2a i)M 2 sinh^)^- 1 ; (60) 

ci = V^{(1 + ^)e- lfcd [cosh(/3 2 ; 2 ) + zX 2 sinh(p^ 2 )] - »(1 - ^)e+ ikd M 2 sinh(/3 2 / 2 )} j-y*-AW , (61) 
^ Pi Pi 

^1 = - ^K 4M [cosh(/3 2 Z 2 ) + iK 2 smh(0 2 l 2 )] - i(l + ^)e +ifed M 2 sinh(/? 2 ? 2 )}J- 1 e( lfc+ ' 3l ) ai ; (62) 

2 pi pi 

ft d6 = e _< * J r " 1 {e"' fcd [cosh(/3 2 Z 2 ) + sinh(/3 2 Z 2 )]Mi sinh(piZi) + e lfed [cosh(/Wi) - iK x sinh(/?iZi)]M 2 sinh(/3 2 Z 2 )}e 2lfeai , 

(63) 

where T = !F[k; /3i,j3 2 ; d; h, l 2 ]. 

VII. APPENDIX B 

In this section, we give a detail derivation of reflective coefficient from the multi-wave interference picture. The 
reflective coefficient as the sum of wave components is 

+oo 

TZ db = ^2n = Ri exp[i(f> rl ] + TfR 2 exp[i(2<£ t i + (j> r2 ) + i2k(d + h)} + T%R%Ri exp{i[2(^ t i + ?V 2 ) + <j>ri] + i2k(2d + h)} 
i=i 

+ T 1 2 i? 2 (i? 2 i? 1 ) 2 cxp{z[(2^ 1 + 4> r2 ) + 2(0 rl + (f> t2 )} + i2k(M + h)} + ... 
1 T?R 2 cx.p[i(2<j> tl + <l> r2 ) + i2k(d+l 1 )] 
" 1 CXp[l<Prl1 + 1 - Rl R 2 eMWn + r2 + 2kd)] 

T^eMWti+klJ-^n+kd}} 

= i?i cxp[i(j) r i\ 



exp(ifcd) — jrrjq exp[— i((f> r i + r2 + kd)] 

„ ! (R 2 /T 2 )T 1 e i ^ exp{# tl - rl + k(2h + d)}} 

= i?i exp [i(j) rl \ , (64 

,F[/c;Pi,p 2 ;d;Zi,Z 2 J 

where in the last step, we used the identities led to eqn. (f2"6")l . Then we can quickly identify 

l + Mfsinh 2 ^) ^ ^ 2 V 7 

i?i exp[# rl ] = ^Mi^ML^ , i? 2 /T 2 = M 2 sinh(/3 2 ; 2 ) (66) 
cosh(piii) + zAi smh(piii) 

by utilizing (|18I19[) and (|20I21[) . Substituting these equations (|65l66p back into we can finally get 

. cosh(,3iii)+'i-Kl sinh(^ih) cosh(/3i!i)-KKi sinh(^iii) C 

d& ~ ~ l Flk-p^fa-d-AM] 

.Mi sinh(/3i/i)J-[fc; ft, 0a; d; l a ] + M 2 sinh(/3 2 / 2 )e ifed 



[ooshOSiii) + itfi smh(^h))T[k; Pi,/3 2 ;d; h,l 2 ] 
.Ah sinh(/3 1 ; 1 )[cosh(^ 2 ; 2 ) + iK 2 smh{f3 2 l 2 )]e- ikd + M 2 sinh( J 9 2 Z 2 )[cosh( y S 1 Z 1 ) - iK 1 sinh(/3 1 / 1 )]e ifed 

F[k;pi,p 2 ;d;li,l 2 ] 

which coincides exactly with (|63l) . 



,(67) 
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(a) Standing Wave in FP Cavity 



Relntivel'rob 




Wiwen limber 



(b) Standing Wave in Vally between Double 
Barriers 



FIG. 7: (a). The relative intensity I(k,x) of a standing wave in FP cavity, with R = 0.80 and d = 2cm. Wave number k is 
plotted with unit 10 6 /m. (b). The relative probability P(k,x) of a standing matter wave in the well of a symmetric double 
barrier, with parameter specified as A — 10.36eV, I = 1.2A and d = 7 A. Wave number is plotted with unit k e as mentioned. 

VIII. APPENDIX C 

In this appendix, we present the analytical expression of phase time of square double-barrier potential. The phase 
shift is defined as 6d bs with 



tan[6>, 



<lh.< 



[M 2 - K 2 cosh(2/3Z)] sm(kd) - K sinh(2/3Z) cos(fcd) 



cosh(2/3Z) cos(fcd) + K sinh(2/3Z) sm(kd) 

So the derivative of 9db s with respect to wave number k is 

d9 dbs <9tan[6> dbs ] 2 T>[k,P,l,d] 

-/[i + tan[Udbs\ \ ~ 



Ok 



Ok 



where S[k,/3,l,d\ = \F[k, (3, 1, d}\ 2 in ^3SJ and 



S[k,/3,l,d\ 



M 



T>[k,P,l,d] = d[M 2 cosh(2/?0 - K 2 ] + — [M smh(2(3l) + K(kl)][l + — (cosh(2/30 - 1)] + 

p I 

2M 2 ( M 2 Kkl 

— — I (cosh(2/?0 - 1)[(1 )sinh(2/3Z) ] cos(2M) + [A'cosh(2/3Z)(cosh(2/3/) - 1) 

P L 2 2 



Thus the analytical expression of tunneling time in the stationary phase approximation is 



(68) 



(69) 



L^]sin(2M)j. 

(70) 



mV\kMd\ 

1~dbs = "T ; J-dbs- 



(71) 



[1] D. Hunger, T. Steinmetz, Y. colombe, C. Deutsch, T.W. Hansch and J. Reichel, New J. Phys. 12, 065038 (2010)[arXiv: 

physics/1005.0067]. 
[2] T. J. Thornton, Rep. Prog. Phys. 57, 311 (1994). 

[3] C.de C. Chamon D.E. Freed, S.A. Kivelson, S.L. Sondhi and X.G. Wen, Phys. Rev. B 55, 2331 (1997) [arXiv: 
|cond-mat/9607195| . 

[4] F.E. Camino, Wei Zhou and V.J. Goldman, Phys. Rev. B 76, 155305 (2007); B.I. Halperin, A. Stern, I. Neder and B. 

Rosenow, Phys. Rev. B 83, 155440 (2011); D.T. McClure, W. Chang, CM. Marcus, L.N. Pfeiffer and K.W. West,Phys. 

Rev. Lett. 108, 256804 (2012). 
[5] Jin-yan Zeng, Quantum Mechanics, (Science Press, 4th edition, 2007), see section 3.3.2. 



14 
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FIG. 8: Wave number is plotted with unit fc e as mentioned above, (a). The 3-Dim graph is specified with barrier height 
A — 10.36eV. (b). The 3- Dim graph is specified with barrier length / = 1.2 A. 
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